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Abstract 

In this work I study the density change along the accretion and wind branches when 
modified Chaplygin gas is taken to be the accreting fluid. I get a high wind density very 
f^ ' near to black hole while a low wind density far from the black hole. In modified Chaplygin 

Cn . gas accretion case a local minima in accretion-density is got which can be detected as the 

r^ ■ CENBOL. 
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The first steps in the field of study about accretion was taken in the five papers of Hoyle, F., Lyttleton, 

R.A. and Bondi, H. [3 [51 131 IH [S] where they have started to analyze the accretion of intersteller fluids 

upon stars in the Newtonian gravity as background. Another particular thing about the Bondi accretion 

was that it was a spherical one. Michel, F. first time has taken the compact object as the Schwarzschild 

singularity and tried to study the accretion procedure as a general relativistic phenomenon. In his study 

the idea of critical point for the first time was noticed and critical conditions were found out. First time 

the approach twoards binary systems, i.e., disc accretion inspite of a spherical one was suggested by 

l^ ■ Shakura, N. I. and Sunyaev, R. A. [7]. Novikov, I. D. and Thorne, K. S. and Page, D. N. [51 [5] enriched 

00 '. the literature by studying the disc accretion in general relativistic perspective. Das, T and Chakraborti, 

^^ ' S. [11] again have shown that as matter is accreting towards a black hole the centrifugal force increases 

^^ I rapidly as compared to the gravitational force. Hence, matter feels increasing centrifugal force around 

t~^ ' a radius little away from the black hole, then slows down and piles up around the radius. This region 

^D I refers to CENBOL (CENtrifugal pressure supported BOundary Layer) acting as the effective boundary 

layer of the black hole system which, like stellar surface, could produce outflowing winds. Therefore, 

due to the presence of CENBOL, a black hole system, like star, easily exhibits wind of matter along 

with inflow. However, once matter crosses the CENBOL, the dominant gravitational force solely controls 

k> ', dynamics and matter falls into the black hole. The nonlinearity arose due to the inclusion of general 

^ • relativity is simplified by deriving the pseudo Newtonian potential. Mukhopadhyay, B. [12] has derived 

5t ! 3' form of pseudopotential for rotating black holes. Using such a potential Mukhopadhyay, B. analysed 

the accretion phenomenon upon black holes where adiabatic gas was the accreting fluid [13). 

Now-a-days there is no dispute that our universe is going through an accelerated expansion phase. 
To explain such an universe we have either to change the left hand side, i.e., the geometry part of Einstein 
equation by modifying Einstein gravity or to change the right hand side's matter part under the Einstein 
gravity. Such exotic matters which can create negative pressure are been differentiated into two types : 
dark energy (DE) and dark matter (DM). Many models are been proposed to describe DE. Chaplygin 
gas model is such a model. Benaoum, H. B. in 2002 has proposed an equation of state (EoS) for Modified 
Chaplygin Gas(modified Chaplygin gas) as 

p^'^p--^ (1) 

As the universe is assumed to be filled up by DE, it is very much obvious that the impact of accretion upon 
black holes would be a matter of interest. In 2004 Babichev, E. et. al. [TS] studided the accretion of dark 
fluid upon a spherically symmetric non rotating black hole. Now besides comparative small black holes 
there exists supermassive galactic black holes. Galaxies centering those black holes are mainly rotating. 
Though the origin of this rotation is quite unknown still. Some assumptions give the responsibility upon 
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primordial turbulance. Consequently, the angular momentum of dark matter halos and eventually the 
rotation of galaxies is thought to be produced by gravitational tidal torque 16 . It has been suggested 
that the halos obtain their spins through the cumulative acquisition of angular momentum from satellite 
accretion. DE that has accreted on a galaxy would be similarly torqued by such tidal interactions. 
Although it is difficult to estimate this effect without invoking specific models, it is reasonable to expect 
that some angular momentum might reside in the DE halos of galaxies. If such rotating DE were to 
accrete on a compact object then it would carry a part of its angular momentum with it, thus leading to 
the formation of a DE accretion disc as opposed to Bondi accretion. 

Biswas, R. et. al. recently [T^ have studided both the spherical and disc accretion upon both rotating 
or nonrotating black holes and analysed the accretion-wind velocity profiles. They have considered steady 
state accretion of inviscid fluid where specific angular momentum is constant throughout the flow and 
got the radial energy momentum conservation equation as 

du 1 dp X^ , ^ 

where all variables are expressed in dimensionless units as follows: 

u — Ur ^ - , X = —, rg = ^^, where Cs^ M and c are the speed of sound, mass of the black hole and the 

speed of light respectively and r, v are dimensionful radial coordinate and radial velocity respectively. 

Here p and p are the dimensionless isotropic pressure and density of the flow, ¥„ — -^ ^, is 

the gravitational force corresponding to the pseudo-Newtonian potential |12| and A is the dimensionless 
speciflc angular momentum (in units of ^^) of the flow, j is the dimensionless specific angular momentum 
of the black hole (Kerr parameter). In the case of Bondi flow A = throughout. 
The equation of continuity is given by 

^ (^«S) = 0, (3) 

where E, the vertically integrated density, is given by |18| E — IcPeh{x), when, Ic = constant (related to 
the equation of state of the accreting fluid) |13| . p^ = density at equatorial plane and h{x) — half-thickness 
of the disc. Assuming the vertical equilibrium, the expression for h(x) can be written as: h{x) ~ Cg 

where c^ = g^ ^ - is the square of sound speed. 
Integrating Q we have 

xi {x^\^ ( nP \(^) 

where is a geometric constant, depending on the exact geometry of the flow. Taking log and differen- 
tiating the equation for sound-speed-square's gradiant becomes 



{l-n)cl + a{n + 1) dcg _ "i I dF 1 du 

{n + 1) Cs (Cg — a) dx 2x 2F dx u dx 

From the EoS the pressure gradiant can be calculated as 

1 dp 1 d , 2\ f a \ d 

p dx n-\- 1 dx \: 

replacing in ([2| and with the help of ([6]) the readial velocity gradiant becomes 

A^ _ p ( \ I /'i Lf^^ E^ 

du im^ ^9\-^) ^ yx Fg dx J {(l-n)c^+a(n+l)} N 

dx ^ u ., ^f" " D' 

u{{l-n)cj+a{n+l)} 



(5) 



Th^{ln(c?-«)}. (6) 



(7) 



Equation ^ has two branches of solution, accretion (inflow) and wind (outflow). As we are considering 
a transonic flow so there will be a point where both N = D = O)^. Using L'Hospital's rule we can obtain 



^Giving the expression for critical point velocity or sound speed immediately as a solution of two different 
algebraic equations. 



the velocity gradiant at that critical point. In the reference jl7| the physically important velocity profiles 
are of adibatic and Chaplygin gas flow of Bondi type, disc flow upon Schwarzschild black hole and Kerr 
type black hole with rotational parameters j = 0.5 and 0.9 respectively. Here I have given these curves 
in the figures 1(a) to 1(h). 

From different literature we can collect the observational data supported values of a and n as 



Table I: 


Permissible 


ranges 


of modified Chaplygin gas parameters in different 


models. 


Source of data 




Type 


of Model 










CDM 


UDME 








a n 


a 


n 


H{z)~z 






0<a< 1.07 < n < 1 


0<a< 1.35 


< n < 1 


Ho — to 
(best-fit) 






a = 0.01 n^O.Ol 


a ^ 0.06 


n^O.ll 



In l(i)-(l) the velocity profile for best fit parameter range from the above table. 

We know for spherical accretion 

M = Anpvr^ (8) 

where M is the change in black hole mass, and the Edington mass accretion rate tells us that 

M is the solar mass= 1.989 x 10^^ gm. 

Now using these two equations if we want to calculate the original fluid density in the accretion flow 
then we must take the values of the units chosen in account. So we get 

GM 
r = X. — TT- ana v — u.c 



So we get the expression of p as 



10^'^c^ 111 

(10) 



Now, as we know 

G = 6.67 X 10-^ gm-^cm^s-^ , c = 3 x lO^^cms-^ 

and letting M — IOMq ,i.e., taking the mass of the black hole as ten times of the solar mass we get, 

p = 1.2214 X 10"^ X At (11) 

UX'^ 

The density-accretion rate equation for disc acretion is given by (|4]) which with the help of (|TT|) 
becomes 

^ 1018 Cs UX^ ^ ' 

Now we have used velocity profiles [T7] of fig 1(a) to 1(1) and plotted p vs x. These are given in 2(a) 
to 2(1). 

If we look at fig 2(a) at first where we have drawn the density profile for spherical accretion of fiuid 
of adiabatic nature upon Schwarzschild black hole. Here the fluid density along accretion branch slowly 
increases and is maximum at black hole event horizon . But in wind branch density increases very slowly 
when radial distance x > Xc, the sonic point distance. But when x < Xc the density increases abruptly 
in the order of . 10^ — 10^. In 2(b) the spherical accretion of modified Chaplygin gas upon Schwarzschild 
black hole has been studied. Here the density in wind branch is steeply increasing as we go near to the 
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Fig. l(a)-l(l) represent the velocity profiles for different paranietrical values. The solid lines 
represent the accretion whereas the dotted lines are for wind. 
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Fig 2(a)-(l) represent the variation of the fluid density in accretion and wind flo"ws as functions 
of radial coordinate corresponding to the velocity profiles sho"wn in (a) Fig 1(a), (b) Fig 1(b), 
(c) Fig 1(c), (d) Fig 1(d), (e) Fig 1(e), (f) Fig 1(f), (g) Fig 1(g), (h) Fig l(hc), (i) Fig l(i), 
(j) Fig l(j), (k) Fig l(k) and (1) Fig 1(1). The solid lines represent the accretion whereas the 
dotted lines are for wind. 



black hole. But the density in accretion branch is bit interesting . At first the density is decreasing as 
we go towards the black hole in the x < Xc region. Then after reaching a local minima again the density 
increases but the rate of increasing is much lower than that of the wind branch. While to analyse 2(a) 
we can tell that the strong gravitational attraction of the black hole makes the fluid more denser as it 
accretes towards the black hole. It is to be noted that the density of wind fluid is much larger than that 
of accreting fluid. As at event horizon the fluid which are to be wind out has the velocity almost equal 
to zero as the strong gravitational force does not allow it to move out. So there the fluid is only gathered 
causing the increment of its density. Besides accreting fluid is been absorbed by the black hole which 
doesn't let the density to be increased. 

On the other hand the case where modified Chaplygin gas is the accreting fluid under the same 
conditions as before we can see the slope of winding fluid is much much higher. We can explain it as we 
go outwards we get a negative pressure so as we go outwards we get a negative prewssure which decreases 
the value of density rapidly. To explain the accreting fluid density we can say that the accreting velocity 
at very far is very low as the negative pressure created by modified Chaplygin gas does not allow it to 
move inwards. So as the fluid is not moving inward or outwardhas a high density there. But as the fluid 
starts to accrete where gravitational pull is stronger than the negative force we will have a lower density 
but when again we are near to black hole again the density increases. 

2(c), 2(e) and 2(g) has almost the same nature as 2 (a) (2(c), 2(e) and 2(g) represent disc flow for 
j=0, 0.5 and 0.9 respectively). Besides the main features of 1(d), 1(f) and l(h)(disc accretion of modified 
Chaplygin gas for j=0, 0.5 and 0.9) are quite similar with 2(b). The only difference to be noted that when 
we go far from the black hole the density is gradually lower for j = case > j = 0.5 case > j — 0.9 case. 
This is very obvious again as increasing black hole spin creates increasing gravitational force which again 
increases the velocity of accretion which does not allow the fluid to gather at a place and force to move 
and the density decreases. Finally, it is to be noted carefully that the modified Chaplygin gas fluid density 
in 2(b),(d),(f),(h) are lower than the corresponding adiabatic cases 2(a),(c),(e) and (g) respectively. This 
is strongly supported by the negative pressure creating nature of modified Chaplygin gas. 

At last we have reached to 2(i) to (1) which are pes x graph corresponding to l(i) to 1(1), i.e., the 
best fit values of a and (3 which are supported by observational data. They all have the basic features as 
2(b). 

Dark energy accretion weakens the accretion and strengthens the wind such that the accretion disc 
may become more and more weak : This was the conclusion of ]17_. In this paper where I am analysing 
the density profiles I can always see that the density along the wind branch is much much higher in 
Chaplygin gas case rather than the adiabatic case. This result resembles with the idea from [17 . But 
the remarkable information in this case is in the accretion curve. Where, accretion branch has a local 
minima showing a minimal value of density near the black hole. It is to be noted thatthis minima is not 
at the critical point. Particularly, I think this local minima signifies the CENBOL. From denominator of 
([7]) it is clear that from D = I do not have Uc — Csc, rather Uc also depends on different parameters like 
a and n etc. According to the literature, beyond CENBOL matter again flow with the supersonic speed. 
For adiabatic cases (1(a), (c), (e) and (g)) it is very clear that identifying the radius where the CENBOL 
is not prominent. But for Chaplygin gas accretion case I can easily identify the radius (let xcenbol) 
where the CENBOL actually takes place and particularly I can particularly say that xqenbol > x^. 
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